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We present thermal measurements of large area (over 1,000 Atm^) monolayer graphene samples at 
cryogenic temperatures to study the electron-phonon thermal conductivity of graphene. By using 
two large samples with areas which differ by a factor of 10, we are able to clearly show the area 
dependence of the electron-phonon cooling. We find that, at temperatures far below the Bloch- 
Gruneisen temperature Tbg, the electron-phonon cooling power is accurately described by the T'* 
temperature dependence predicted for clean samples. Using this model, we are able to extract a 
value for the electron-phonon coupling constant as a function of gate voltage, and the graphene 
electron-lattice deformation potential. We also present results for thermal conductance at higher 
temperatures, above Tbg/4, for which the clean limit no longer applies. In this regime we find a 
cooling power which is accurately described qualitatively, but not quantitatively, by a model which 
predicts the emission of very high energy phonons through a disorder-assisted mechanism. 


I. INTRODUCTION 

The potential applications for graphene as a highly sen¬ 
sitive photon detector have driven substantial interest in 
determining the thermal conductance of graphene at low 
temperatures [1, 2]. In addition to the scientihc value 
of a greater understanding of electron-phonon coupling 
in graphene, knowledge of this important physical pro¬ 
cess is critical to determining the theoretical performance 
of highly sensitive graphene-based photon detectors. If 
the thermal conductance of an ultra-sensitive graphene- 
based photon detector is too large, the detector will cool 
off too quickly to allow for accurate photon detection 
[3, 4]. To date, thermal measurements of graphene have 
used Johnson noise thermometry [1, 5-7], the temper¬ 
ature dependent resistance of a superconducting tunnel 
barrier [8] or supercurrent hysteresis [9] to measure the 
thermal conductance of graphene as a function of dissi¬ 
pated power, which determines the electron temperature. 
However, due to wide variability between samples, even 
within a single study, it has been difficult to form gen¬ 
eral conclusions about the phonon cooling pathway in 
graphene. In addition, the measurements are challeng¬ 
ing. Typically the sample resistance is large, the contacts 
can add extra resistance and provide an additional cool¬ 
ing pathway, and the signals in these measurements are 
small. We discuss these issues below in detail. 

In metallic thin films (film thickness ^ 10 nm to 1 /rm), 
the electron-phonon cooling power typically takes the 
form P = VTi (T^ — Tq ), where V is the device volume, 
S is a coefficient that describes the electron-phonon cou¬ 
pling strength, and J is a parameter that varies from 
4 — 6, dependent upon the amount of disorder in the sys¬ 
tem [10]. Similarly, in monolayer graphene the electron- 
phonon cooling power and thermal conductance, Gep, de¬ 
pend on the level of disorder, screening, and the temper¬ 
ature of the electron system. We discuss the phonon 


cooling in the applicable theories, and then present our 
experiments on two monolayer graphene samples. 

By accounting for the microwave losses associated with 
substrate electrons, we are able to accurately measure 
the cooling of the samples. We find that, at low temper¬ 
atures, the cooling in both samples can be described by 
a combination of hot electron out-diffusion and phonon 
emission. By tuning the carrier density of the samples 
and measuring devices with areas that differ by an or¬ 
der of magnitude, we are able to definitively identify 
the phonon cooling pathway. We find Ggp to be con¬ 
sistent with the behavior expected for a pure monolayer 
graphene sheet with little disorder. 


II. THEORY OF PHONON COOLING 
A. Clean limit 

The cooling power in the limit of large electron- 
impurity mean free path ^mfp (taken to be inhnite) fol¬ 
lows a power law form 

P = - T^) (1) 

and is due to the emission of longitudinal acoustic (LA) 
phonons into the graphene lattice [11]. In Eq. 1, T and 
To are the temperatures of the electron system and the 
lattice, respectively. The forms of the coefficient E and 
exponent <5 depend on the temperature of the electron 
system, with a crossover temperature of order the Bloch- 
Gruneisen temperature Tbg = 2sy^^/fcB [12], where 
s = 2 X 10^ m/s is the speed of sound in graphene, n 
is the carrier density in graphene, and Icb is Boltzmann’s 
constant. For n = 10^^ cm“^, as is typical for samples 
on Si02, Tbg = 54 K. At temperatures T ^ Tbg, 

P = AEi(t4-T4), 


( 2 ) 
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where 

^ Ibpufi^VpS ^' 

Here, pu is the mass density of graphene, vp = 10® m/s 
is the Fermi velocity in graphene, Ep = hyp^/im is 
the Fermi energy, and D is the deformation potential of 
graphene. The deformation potential is a measure of the 
strength of electron-phonon coupling and has been stud¬ 
ied theoretically [13-15] and experimentally [5, 6, 12, 16]. 
Calculations find D ranges from 2 to 70 eV, with theoret¬ 
ical predictions ranging from approximately 5 to 13 eV. 

For T ^ Tbg, the cooling power for typical devices, 
with Ep ^ fceT, is given by 

P = 5iH(r-To), 

where 

^ D^E^kn 

2TTpMh^V% ' 

For very low carrier densities or operation at high T, 
where Ep < k^T, 


(4) 

(5) 
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^ 30/^M^®^'F ' 


( 6 ) 


This is not a regime we access in our measurements, as 
Ep is approximately equal to k-^T at T = 1000 K at 
a carrier density of n = 10^^ cm“^. On the disordered 
Si 02 substrates which are used here, much lower carrier 
densities are difficult to achieve [17] and temperatures 
above 1000 K are inaccessible. 


B. Disorder-assisted scattering 


The introduction of disorder (by means of a finite mean 
free path) has ramifications at both low temperatures and 
high temperatures. Theoretical 

[18] find a different form of electron-phonon cooling 
power than given by Eq. 2 at low temperatures. It is 
predicted that at temperatures below a crossover tem¬ 
perature Tx, the electron-phonon coupling is enhanced 
and is larger than the value given by Eq. 2: 

p = HE2(r®-r3), (7) 

where 

_ 2C(3)P^lPF|fc| . . 

and C{n) is the Reimann zeta function; C(3) « 1.2. The 
crossover temperature Px is the temperature for which 
the expressions in Eqs. 2 and 7 are equal. This temper¬ 
ature is given by 


Px 


30fcC(3) 

TT^/cBt^mfp 


(9) 


Using a Drude model, we find a lower bound of £mfp = 
50 nm for our samples, which sets an upper bound of 
Px « I K. 

At higher temperatures, above PbGj a different form of 
electron-phonon scattering is predicted [19] in which dis¬ 
order in the graphene allows large momentum phonons 
to be emitted from the electron system. These large 
momentum scattering events dissipate energy of or¬ 
der fcuP and are coined “supercollisions.” In typi¬ 
cal momentum-conserving scattering, the momentum of 
scattered phonons is constrained by the Fermi surface 
of the graphene (a circle in fc-space of radius \Ep\/hvp). 
However, supercollisions allow phonons with much larger 
momenta to be emitted, with an accompanying recoil 
phonon which allows the net energy exchange to be 
momentum-conserving. The result is that the cooling 
power due to supercollisions for P > Pbg is predicted 
to be exactly half the cooling power predicted for low 
temperature disorder-assisted scattering: 

P=^A(p3-P3). (10) 

In order to predict the electron-phonon thermal con¬ 
ductance in graphene, it is critical to clearly state the 
assumptions made about the presence of disorder and 
the temperature regime in which the device is operating. 
Moreover, measurement of Gep is nontrivial. Though 
several groups have made measurements of graphene’s 
thermal properties at cryogenic temperatures [1, 6, 9], 
there is considerable disagreement between experiments 
about the functional form of Gep as well as the mag¬ 
nitude of D [11]. In our measurements, we study very 
large area graphene samples to emphasize Gep, which 
is proportional to the graphene area. This enables 
us to clearly separate the electron-phonon cooling in 
our samples from other cooling pathways (notably, hot- 
electron out-diffusion). In addition, by looking at similar 
graphene samples with substantially different area, we 
are able to accurately probe the area dependence of the 
electron-phonon cooling channel. 


III. DEVICE FABRICATION 

The samples used in our measurements are prepared 
from CVD-grown graphene purchased from ACS Mate¬ 
rial [20]. The graphene is grown using copper foil as 
a catalyst. Prior to purchase, the graphene was trans¬ 
ferred to an oxidized doped silicon substrate by coat¬ 
ing the graphene sample with PMMA and etching away 
the copper foil. It has been shown that this process can 
leave PMMA residue on the surface of the graphene [21]. 
However, we found that the ability to lithograph large 
graphene areas to emphasize Gep and reduce contact re¬ 
sistance outweighed the negative aspects of having pos¬ 
sible contaminants. The doped Si substrate with a room 
temperature resistivity of p = 1 il • cm allows for the car¬ 
rier density of the graphene to be controlled in situ by 
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Sample 

L (/im 

l) W (/tm) 

R (Q) 

Rn m 

D (eV) 

Gl 

10 

100 

90 

900 

12.0 

G2 

100 

100 

840 

840 

10.5 


TABLE I. Comparison between the shorter and longer 
graphene samples (G1 and G2, respectively) at a gate voltage 
of 14 « —32 V. The values of D were obtained by individ¬ 
ual fitting the thermal conductance as a function of T using 
Eq. 12 for T < Tbg/4^ at each gate voltage, Fig. 6 below. 

using an electrostatic gate voltage to capacitively induce 
either electrons or holes. However, the doped substrate 
can complicate the measurement of device temperature 
(discussed in Section IV). 

After purchase, the samples are fabricated in a multi- 
step process. Using electron-beam lithography, we define 
areas of the graphene sheet which we then remove with 
an oxygen etch. What remains are graphene sheets of 
width 100 /im and lengths ranging from 10 to 200 fim. 
Then, we again use the electron beam patterning, this 
time to define the contact and lead structure. We de¬ 
posit a Pd/Al bilayer with thicknesses 5 nm/50 nm to 
contact the graphene. Palladium is used to form Ohmic 
contact to the graphene [22, 23] and aluminum was used 
to realize superconducting contacts with a T(. « 1.2 K. 
Superconducting contacts were desired to suppress the 
out-diffusion of hot electrons, another potential source of 
thermal conductance to the bath. 


both with channel width W = 100 /rm. The first has 
length L = 10 ^m between contacts (sample Gl) and 
the other has a channel length of 100 /im (sample G2). 
Measuring these two devices in a single cool down of 
the refrigerator allowed us to establish how much con¬ 
tact resistance is present and study the area-dependence 
of thermal conductance. In Fig. 2 we plot the resistance 
per square (!?□) of both devices as a function of gate volt¬ 
age. We define Ra = R{W/L) where R is the measured 
resistance. The measured resistance can have a contri¬ 
bution due to the series contact resistance if present, but 
none due to the Al/Pd leads for T < Tc. At negligible 
contact resistance, R should scale with length for fixed 
W. Indeed, we find that R\j is approximately the same 
for the two samples. At gate voltages jVgj > 10 V, far 
from the charge neutrality point (CNP), we calculate a 
mobility of approximately 3500 cm/s^ for both electrons 
and holes, which is consistent with high quality GVD- 
grown graphene on Si02 [24]. We find that the resis¬ 
tance is relatively insensitive to bias current over several 
orders of magnitude of current. The comparison of R\j 
at Vg = —32 Y for samples Gl and sample G2 indicates 
that the contact resistance is small compared to the total 
device resistance. In addition, analysis of the conductivi¬ 
ties as a function of gate voltage [25] yields residual root- 
mean-square carrier densities of n « 2.5 x 10^^ cm“^ and 
3 X 10^^ cm“^ at the resistance peak for samples Gl and 
G2, respectively. This is due to carrier puddling arising 
from external electric fields [26]. 


IV. MEASUREMENTS 

The measurements for this study were performed in 
an Oxford Triton 200 cryogen-free dilution refrigerator. 
With no cabling, this refrigerator can achieve a base tem¬ 
perature of 8 mK. 
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FIG. 1. Scanning electron microscope image of a typical 
graphene sample. The metallic leads are false colored in yel¬ 
low. The dark regions are areas of graphene which were not 
etched away and the remaining light area is the Si 02 sub¬ 
strate. Note that the channel length L for the sample in the 
image is 20 /tm. For the samples measured, L is 10 or 100 /tm. 
For both the sample shown in this figure and those measured, 
the channel width W is 100 /tm. 

We present the thermal measurements of the two sam- 
ples(typical geometry shown in Fig. 1) with device prop¬ 
erties given in Table I from the same graphene growth. 


FIG. 2. Resistance per square, Ra, as a function of gate 
voltage for sample Gl (L = 10 /tm) and sample G2 (L = 
100 /tm) at a bias current of / = 1 /tA. The gate voltages 
of each curve are offset so that the resistance peak occurs at 
14 = 0. The offsets are —1.5 V and —8 V for samples Gl and 
G2, respectively. Inset: Resistance per square as a function of 
bias current at I 4 « —32 V. Data were taken at To < 100 mK. 

To study the thermal properties of the graphene, a 
constant power is applied to the electron system and the 
emitted Johnson noise of the system is measured to deter¬ 
mine the resulting change in electron temperature. These 
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FIG. 3. Apparatus for Johnson noise measurements. The 
metallic leads contacting the device are represented by the 
hatched regions. The low-pass and band-pass filters are indi¬ 
cated by LP and BP (with center frequency / = 1.3 GHz), 
respectively. The bias-T at base is represented by a capaci¬ 
tor and inductor. The rf signal from the spectrum analyzer 
is attenuated by 40 dB before coupling to the device and is 
necessary for the reflectometry measurement described in the 
text. 


power is also plotted. Here we consider the impedance 
mismatch between the 50 H microwave line and an equiv¬ 
alent resistor with the dc resistance of the graphene sam¬ 
ple. Although the rf impedance at 1 GHz of the graphene 
itself is approximately equal to its dc resistance [27], 
the data deviates substantially from the calculation seen 
in Fig. 4a, where the dc resistance is used to calculate 
the circuit rf impedance. We believe the discrepancy is 
not due to the graphene, but instead arises from charge 
carriers in the doped Si substrate capacitively coupling 
through the large contact pads. Using the very simple 
circuit model shown in Fig. 4c where the carriers in the 
substrate provide a parallel resistance, we can approxi¬ 
mately replicate the observed microwave behavior, with 
the values of the lumped circuit elements given in the 
figure caption. 


steady state measurements allow us to probe the cool¬ 
ing pathways of the device. We use a Yokagawa 7651 
voltage source and a large bias resistor to apply a dc 
or low frequency on/off current of amplitude I to the 
graphene (see Fig. 3). The current heats the electrons 
in the graphene with power P = I^R. This heating 
power raises the steady state electron temperature of the 
graphene above the stage temperature Tq. This change in 
electron temperature results in a change in emitted John¬ 
son noise power into a matched load equal to k^BAT, 
where AT = T—Tq. This result assumes that the emitted 
noise is measured at a sufficiently low frequency / so that 
hf ^ k^T, as is the case in our study with / = 1.3 GHz, 
and that the electron temperature is constant across the 
device. In the case of finite thermal conductance from 
charges diffusing out the leads, the electron temperature 
as a function of position T(x) is not constant, so the 
measurement of Johnson noise probes the average elec¬ 
tron temperature of the graphene T: 

L 

T = — y T{x)dx (11) 

0 


and AT = f — Tq. 

The Johnson noise signal is rectified to produce a dc 
voltage using a zero bias Schottky diode. The change in 
diode voltage is given by AVdiode = kAT, where k is a 
coupling constant representing the amplification of the 
50 Q microwave output system. In order to accurately 
measure the average temperature of the electron system, 
it is thus critical to determine the value of k. This is one 
of the important calibration procedures necessitated by 
the use of a commercially-supplied doped Si substrate, 
which is weakly electrically conducting. 

To better understand the microwave coupling to the 
device, we performed reflectometry measurements using 
a spectrum analyzer. In Fig. 4a the normalized reflected 
power as a function of gate voltage is plotted for sample 
Gl. For comparison, the expected normalized reflected 




Gate Voltage (V) Average Electron Temperature (K) 



FIG. 4. Gharacterization and calibration of Sample Gl. (a) 
Reflected power from Sample Gl as a function of gate voltage. 
The DG model gives the reflected power that would be ex¬ 
pected from calculating the reflection coefficient using the de¬ 
vice resistance. The LCR model instead calculates reflection 
using the circuit model of (c), with L = 2.5 nH, Gsub = 10 pF, 
Rsuh = 290 H, and Rg equal to the internal graphene resis¬ 
tance. The capacitance was calculated from the area of the 
lead structure, the inductance was estimated from wirebond 
length, and the substrate resistance was adjusted to provide 
a good fit to the data. At V), = —32 V, the circuit model 
yield a power reflection of F^ = 0.15, while the DC model 
predicts F^ = 0.09. Thus, the device is well matched when 
gated far from the CNP. (b) Calculated thermal conductance 
as a function of average electron temperature for several base 
temperatures. To determine the coupling, k is chosen (ex¬ 
plained in the text) to align the thermal conductance curves 
(so that G{T) is independent of Tq). 
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The coupling of substrate charges to the Al/Pd leads 
means that the carriers in the substrate also emit John¬ 
son noise, which is measured at the output. As a re¬ 
sult, a calibration based on the diode output voltage as 
a function of stage temperature would yield a slope that 
is larger than k. Instead, we determine the coupling by 
measuring the diode voltage as a function of DC input 
power at different stage temperatures and requiring that 
the G{T) = dP/dT be independent of the stage tem¬ 
perature at which the measurement was performed (see 
Fig. 4b). It is important to note that this method is only 
exactly correct in the case of uniform electron tempera¬ 
ture across the graphene device. In our system, where 
we observe a non-uniform temperature distribution due 
to carrier out-diffusion, we find from simulations that this 
calibration technique can understate the value of k by 10 
to 15%. We account for this effect in the data we present 
in Figs. 5,6, and 7. 


A. Cooling model 

In analyzing the thermal properties of the samples, the 
heat diffusion equation is used to model to cooling power 
of the system: 

Cr = p.p(x)-A(„M^), (12) 

where r is the resistance per unit length, Pep{x) is 
the electron-phonon cooling power per unit length, and 
g{x) = CT/R [2 is the thermal conductance from carrier 
diffusion. The appropriate form of Pepix) depends both 
on the electron temperature and on the level of disorder 
in the system (i.e., £mfp)- For example, in the low tem¬ 
perature, clean limit, Pep{x) = WSi {T{x)'^ — Tp). After 
solving for T(x) at a given current, the average electron 
temperature of the system T, can be found by integrating 
over the entire length, as per Eq. 11. 

This average temperature is what is measured using 
the Johnson noise method. Note that this calculation is 
correct only if the resistance is approximately tempera¬ 
ture independent, which is very close to correct for our 
samples (Fig. 2). 

The boundary condition used to solve Eq. 12 is that 
T{x) = Tq at both ends of the device, as there was strong 
evidence of diffusion cooling of hot electrons in sample 
Gl. Because sample G2 has a resistance and area that 
is approximately 10 times larger than that of sample Gl, 
diffusion cooling is not evident in our experiments. The 
measurements of T were performed over a large range 
of excitation powers (/ = 0.1 — 1000 pA), so we access 
electron temperatures above and below Tbg in our mea¬ 
surement. We model these two regimes separately, first 
focusing on low temperatures, below Tbg/4. 



FIG. 5. Thermal conductance, G = dP/dT of sample Gl. 
The solid line is a calculation of thermal conductance for 
T < Tbg/J using the device parameters given in Table II and 
Eqs. 2 and 12. The same calculation is shown by a dashed 
line for T > Tbg/4. Inset: calculation of Ei using n calcu¬ 
lated from Vg and measured values of D (given in Table II). 
The dashed line shows the value Ei as a function of n for 
D = 11.5 eV, demonstrating the expected Ei oc ydi depen¬ 
dence. The uncertainty at each of the data points for Ei arises 
from the calibration procedure at each gate voltage. There 
is an additional overall uncertainty of scale for Ei (of ±20%) 
due to uncertainty in the stage thermometer calibration. 


(V) 

n (10^^ cm' 

■^) Tbg (K) 

R (fi) 

D (eV) 

-33.5 

24.1 

84 

90 

12.0 

-25.5 

18.4 

73 

110 

11.0 

-17.5 

12.6 

60 

150 

11.0 

-9.5 

6.8 

45 

250 

11.5 

-3.5 

2.5 

27 

580 

11.5 


TABLE II. Device properties of sample Gl {L — 10 pm). The 
values of D were obtained by individually fitting the thermal 
conductance as a function of T using Eq. 12 for T < Tbg/J 
at each gate voltage, Fig. 5. 


B. Low temperature cooling 

In Fig. 5, we plot the thermal conductance of sam¬ 
ple Gl for T < 20 K as a function of the average elec¬ 
tron temperature T for multiple gate voltages (in the hole 
doped region). The solid lines are plotted for T < Tbg/ 4 
for each density. These plots are a calculation of the effec¬ 
tive thermal conductance, G = dP/dT, using the temper¬ 
ature distribution determined by Eq. 12. At low temper¬ 
atures, for each gate voltage, the thermal conductance of 
this sample is approximately equal to the value predicted 
for the out-diffusion of hot electrons to the contacts a 
one-dimensional wire: Gdis ~ liCT/R, where C is the 
Lorentz number [28, 29]. In our measurements, we find 
£ = 1.18 X £theory, where £theory = 2.45 X 10“® WGK“^. 

At temperatures above a few Kelvin, the cooling is in- 
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creasingly dominated by the emission of phonons, and we 
use the clean limit form of Pep to calculate the electron 
temperature. This determines Si. The only free pa¬ 
rameter in determining the thermal conductance is the 
deformation potential which appears in Eq. 2. For 
each gate voltage, the value of D was independently de¬ 
termined from Si by fitting to the relevant data points, 
and the plotted lines represent the best fit. The values 
of D found range from 11 to 12±2.0 eV and are pre¬ 
sented along with other relevant physical parameters in 
Table II. With these values of D, we have calculated S 
as a function of n and have plotted these results in the 
inset of Fig. 5. The dashed line in the inset represents the 
anticipated n dependence of Si oc ^/n for a deformation 
potential of D = 11.5 eV. 

We also tested the disorder-assisted scattering model, 
with Pep given by Eq. 7, in solving Eq. 12. We found that 
this model did not agree with the data. We conclude 
that the data do not exhibit an electron-phonon cooling 
power with a dependence at temperatures below 20 K. 
This conclusion is consistent with measurements of device 
resistance, which indicate a disorder-limited mean free 
path of over 50 nm. For this value of £mfp, the crossover 
temperature T^, below which disorder-assisted scatter¬ 
ing would be expected to play a role, is Tx « 1 K. For 
all gate voltages measured, the thermal conductance at 
these lowest temperatures (T < 1) is dominated by elec¬ 
tron out-diffusion. Thus, we could not test the disordered 
limit in the present experiments 

The clean limit calculation ceases to agree with the 
measured data at large temperatures, and the temper¬ 
ature at which this occurs is dependent upon the gate 
voltage. This is to be expected, as Tbg cic i/n so that 
Tbg/ 4 « 11 K at Vg = —9.5, for example, whereas 
T’bg/ 4 « 21 K for Vg = —33.5. Thus, the tempera¬ 
ture range in which the low-temperature clean limit ac¬ 
curately describes Pep is smaller for lower carrier densi¬ 
ties. 

We tested the predicted area dependence of Eq. 2 by 
comparing the thermal conductance of the shorter and 
longer samples (G1 and G2, respectively). In Fig. 6, 
we plot the thermal conductance at high density, Vg « 
—32 V, for both sample G1 and sample G2. At low tem¬ 
peratures, the total thermal conductance of sample G2 
is less than that of sample G1 due to the approximately 
10 times larger internal resistance of sample G2 (see Ta¬ 
ble I). The increased length of G2 suppresses Gdis by a 
factor of 10, resulting in a much lower total thermal con¬ 
ductance below 1 K. However, at high temperatures, the 
thermal conductance of sample G2 is approximately 10 
times larger than that of sample G1 as expected. The fits 
for samples G1 and G2 yield D = 12 eV and D — 10.5 eV, 
respectively, in good agreement and consistent within the 
measurement uncertainty. 



Average Electron Temperature (K) 

FIG. 6. Plot of thermal conductance for samples G1 and G2. 
The solid lines are calculations of thermal conductance using 
Eq. 12 in the low temperature clean limit, with D = 12 eV 
and D = 10.5 eV for samples G1 and G2, respectively. The 
dashed vertical line indicates Tbg/4 « 21 K. The temperature 
ranges over which either Gdiff or Gep is the dominant source 
of thermal conductance in sample G1 are indicated by the 
arrows in the figure. For sample G2, Gep S> Gdift for all 
temperatures measured. 

C. High temperature cooling 

At high bias currents, the electron temperature can 
exceed Tbg/4 so that the electron-phonon cooling of the 
graphene can no longer be described by its low temper¬ 
ature limit. In Fig. 7 we plot the thermal conductance 
of sample G1 at many gate voltages for T > 10 K. Over 
a temperature range of approximately 15 K to 50 K, the 
thermal conductance follows an approximately T^ power 
law, which suggests that supercollisions may contribute 
to Gep. As a result, we used a linear sum of Eqs. 4 and 10 
to model the electron phonon cooling power, using the 
values for D that were extracted at lower temperatures. 
In these calculations, the only free parameter is the mean 
free path due to short range disorder potentials. 

We find that the best fit to the supercollision model is 
obtained with £mfp = 3 ± 1 nm. This value for electron- 
impurity mean free path is quite low and is much lower 
than would be expected from the lower bound obtained 
by using the Drude model with the device resistance. 
Moreover, such a small value for would suggest that 
at lower temperatures, disorder-assisted scattering would 
have a substantial contribution to the thermal conduc¬ 
tance. Such behavior is not observed. Other groups us¬ 
ing supercollisions to model their data have found very 
large values for the deformation potential {D > 50 eV) 
by fixing to an estimated value [6, 16]. Because the 
coefficient E 2 in Eq. 10 is proportional to T^/^mfp, both 
our results and those of previous studies find that at¬ 
tempting to explain the thermal conductance data with 
the supercollision model results in extracted physical pa- 
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Study 

Regime 

Substrate A (pm^) 

R (kH) 

D (eV) 

Betz et al. [6] 

Supercollisions 

BN 

6.2 

1.5 

70 

Betz et al. [7] 

Clean limit 

BN 

6,13 

1,3 

4,2 

Fong et al. [1] 

Clean limit 

SiOa/Si 

102 

30 

33 

Fong et al. [5] 

Low-T disorder 

SiOa/Si 

25,55 

1.5,5 

51,23 

Theory [13-15] 

- 

- 

- 

- 

5-13 

This study 

Clean limit 

SiOa/Si 

10 "*, 10® 

0.09,0.84 

12,11 


TABLE III. Summary of experimental results which studied the low-temperature electron-phonon coupling in graphene using 
Johnson noise thermometry. Individual samples within a study are separated by commas. 


rameters D and/or ^nifp that are not within the expected 
range for those devices. 

At T > 50 K, the thermal conductance appears to 
increase faster as a function of T than would be expected 
for a dependence. Also, the curves asymptote to a 
value which is independent of gate voltage. This behavior 
is not anticipated from models of LA phonon scattering 
within the graphene, so we believe that this is evidence 
of the emission of optical phonons directly into the Si 02 
substrate [15, 30, 31]. At these high energies, the large 
potential energy drop over a typical electron mean free 
path might exceed the energy needed to emit a surface 
polar photon. 



Average Electron Temperature (K) 


FIG. 7. Measured thermal conductance, G = dP/dT, of 
sample Gl. The solid line is a calculation of thermal conduc¬ 
tance for T > Tbg/J using the device parameters given in 
Table II, ^mfp = 3 nm, and Eqs. 4, 10, and 12 This value of 
^mfp is much smaller than expected, as discussed in the text. 


V. DISCUSSION 

The Johnson noise emission measurements reported in 
the previous section probed the cooling processes of two 
graphene samples with areas which differed by a factor 
of 10. This made it possible to observe the dependence 
of cooling power on the device area which is consistent 
with theory and, in the case of the longer device, remove 


the effects of electron out-diffusion. The measurements 
of both the 10 /rm long and of the 100 /rm long sam¬ 
ple found similar results for the deformation potential 
D, with measured values ranging from 10.5 — 12 eV. Due 
to systematic uncertainties (primarily arising from the 
calibration and thermometry), there is an overall uncer¬ 
tainty in D of approximately ±20%. As a result, the 
bound that can reasonably be placed on the deforma¬ 
tion potential from these measurements is 8.5 — 13.5 eV 
(see Table II). In Table III, we compare this result to 
the extracted values of D obtained from other studies, 
which show wide variation from sample to sample and 
little agreement with the theoretical predictions. 

At low temperatures, the electron-phonon cooling is 
consistent with the clean limit for low electron temper¬ 
atures (T < 20 K for n = 2.4 x 10^^/cm^), for which 
G (xT^ is predicted. We did not observe any behavior 
consistent with low temperature disorder-assisted scat¬ 
tering [18] in the shorter samples, but the cooling path¬ 
way of hot electrons diffusing out the leads might obscure 
a deviation from the clean limit in the shorter sample, 
Gl, at low temperatures if one were present. For the 
much longer sample, G2, a small deviation from G oc 
was observed at very low temperatures for all gate volt¬ 
ages (see Fig. 6). However, there is greater relative uncer¬ 
tainty in the graphene electron temperature below 1 K, 
which precludes drawing any quantitative conclusions for 
the electron-phonon cooling processes at these very low 
temperatures. 

A interesting future study would be thermal measure¬ 
ments of a graphene sample with low electrical resis¬ 
tance, which was contacted with high Tc metallic con¬ 
tacts (Tc > 9 K) that were able to confine hot electrons 
and suppress Gdiff. A low sample resistance would be 
desirable to achieve good microwave coupling to the out¬ 
put lines and to signihcantly mitigate the effects of free 
carriers in the substrate. Clearly, an insulating substrate 
would be desirable for these measurements, with a sep¬ 
arate gate electrode which couples to the sample at dc, 
but not at microwave frequencies. In addition, if the sam¬ 
ple were suspended [25, 32], thermal conductance mea¬ 
surements at high dc bias would provide a useful test of 
whether optical surface phonons can explain the sharp 
upturn in thermal conductance observed for T > 50 K. 
However, care must be taken to ensure that suspending 
the device does not generate a bottleneck for removal of 
excess phonons from the graphene lattice as this could 













complicate the interpretation of the data. 
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